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Abstract. It is shown that to leading order the collisional energy loss of an energetic parton in the hot 
quark gluon plasma reads dE/dx ~ a(mjj)T 2 y where the scale of the coupling is determined by the 
(parametrically soft) Debye screening mass. Compared to previous expressions derived by Bjorken and 
other authors, dE B /dx ~ a 2 T 2 ln(ET /m 2 D ), the rectified result takes into account the running of the 
coupling, as dictated by quantum corrections beyond tree level. As one significant consequence, due to 
asymptotic freedom, the QCD collisional energy loss becomes independent of the jet energy in the limit 
E ^> T. It is advocated that this resummation improved perturbative result might be useful to (re-)estimate 
the collisional energy loss for temperatures relevant in heavy ion phenomenology. 

PACS. 12.38.Mh Quark-gluon plasma, 25.75.-q Relativistic heavy-ion collisions 



1 Introduction 

The quenching of jets observed at RHIC [1] has become 
one of the key arguments for the creation of a 'new state 
of matter' in heavy ion collisions. Being, inter alia, at- 
tributable to the energy loss of energetic partons ('jets'), 
it indeed allows to infer specific properties of the traversed 
matter. In a quark gluon plasma (QGP) there are two dif- 
ferent loss mechanisms: elastic collisions with deconfmed 
partons [2,3], or induced gluon radiation [4-7]. Presuming 
a dominance of the second effect, experimental findings 
have often been interpreted in terms of a purely radiative 
energy loss. However, the data- adjusted model parameters 
(either q or dN g /dy, depending on the approach) seem de- 
batable because they are considerably larger than theoret- 
ically expected or even in conflict with a strong constraint 
from dS/dy, see e.g. [8]. 

An obvious remedy for this situation is a sizable col- 
lisional contribution to the energy loss, as considered re- 
cently in [9,10]. These considerations as well as all pre- 
vious studies were based on formulae of the generic form 
dE B /dx ~ a 2 T 2 \n(ET/m 2 D ) due to Bjorken, which can 
be considered an adaption of the (QED) Bethe-Bloch for- 
mula. Evidently, such estimates depend crucially on the 
size of the coupling a. For practical estimates it has been 
common practice to choose a number for a by evaluating 
the running coupling at a 'typical' thermal momentum 
scale, often of the order of 2ttT. This ad hoc 'prescription' 
was scrutinized only recently [11]. 

In fact, Bjorken's formula is based on a tree level scat- 
tering amplitude - hence the parameter a in dE B /dx is to 
be considered a fixed number. This level of approximation 
is acceptable in the QED analog, where the strength of 
the electromagnetic coupling is virtually constant over a 
huge range of momenta. The strong coupling, on the other 



hand, varies considerably in a momentum range which can 
be probed by the energy loss in heavy ion collisions. Thus, 
in QCD, quantum corrections to the tree level scattering 
amplitude, which induce the running of the coupling, have 
to be taken into account on conceptual grounds. 

Before discussing the implications of this statement for 
the collisional energy loss, I will briefly review the deriva- 
tion of Bjorken's formula. 



2 Bjorken's formula 

Bjorken [2] considered the propagation of a jet j through a 
static QGP at a temperature T 3> A, where the coupling 
is small and a perturbative approach is appropriate. Its 
mean energy loss per length can be calculated from the 
rate of binary collisions with partons of the medium, as 
determined by the flux and the cross section, 
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Here p s 
and d a 



d s n s is the density of scatterers s, with d g = 16 



* q — I2nj being the gluon and quark degeneracies 
for n/ light flavors, and n s — > n±(k) = [exp(fc/T) ± 1]~ 
in the ideal gas approximation. Furthermore, denotes 
a dimensionless flux factor, t the 4-momentum transfer 
squared, and lo = E — E' the energy difference of the in- 
coming and outgoing jet. Scatterings with small momen- 
tum transfer dominate, thus the cross section can be ap- 
proximated by the i-channel contribution, 
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with C qq — |, C qg = 1, and C' gg = |. For E and £" 
much larger than the typical momentum of the thermal 
scatterers, k <~ T, the relation of i to the angle between 
the jet and the scatterer simplifies to 



t = -2(1 - cos0)kuj, 



(3) 



and the flux factor can be approximated by <P = 1 — cos 9. 
At this point, Bjorken integrated in Eq. (1) 
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imposing both an IR and UV regularization. The soft cut- 
off is related to the Debye mass, \t 2 \ = n 2 ~ m 2 D ~ aT 2 , 
describing the screening of the exchanged gluon by the 
medium. The upper bound of \t\ was reasoned to be given 
by the maximum energy transfer: very hard transfers, say 
lo w E, do not contribute to the energy loss; in this case 
the energy is collinearly relocated to the scatterer. As- 
suming LOmax — E/2 implies t\ = —(1 — cos9)kE, thus 
dEf/dx = ira 2 J2 S Qs J k3 k^ps In((l - cos 0)kE/ f i 2 ) . Re- 
placing now, somewhat pragmatically, the logarithm by 
ln(2(fc) E j ' \x 2 ) and setting (k) — ► 2T, Bjorken arrived at 
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which differs for quarks and gluons only by the prefactor. 

Various refinements of Bjorkcn's illustrative calcula- 
tion, as e.g. [3], aim at the precise determination of the 
cut-offs. The generic form of these results agrees with 
Eq. (5), with fi — ► mjj in the logarithm and the factor 
4 replaced by some function of n/. 



3 QCD collisional loss with running coupling 

As already motivated in the Introduction, the running of 
the coupling should be regarded when its strength varies 
notably in the range of probed momenta. In the present 
case of a t-channel dominated observable, the relevant 
scale is the virtuality of the exchanged gluon, t E [t\,t2\- 
For a representative situation (T = 0.3 GeV, E = 10 GeV), 
one roughly finds |*i j ~ 10 GeV 2 and \t 2 \ ~ lGeV 2 , thus 
for phenomenological estimates it is conspicuously essen- 
tial to take into account the running of the coupling. 

Formally, the running coupling arises in the course of 
renormalization of the (divergent) loop corrections to the 
tree level scattering amplitude. Referring to [11] for some 
details of the elementary considerations, the resulting re- 
summation improved cross section is basically the expres- 
sion (2) with the (constant) parameter a replaced by 
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where A is the QCD parameter, and &o = (11— 3 n/)/(47r). 
Moreover, an infrared cut-off fi 2 <~ m 2 D (stemming from 



the thermal contributions to the gluon selfenergy) is to be 
understood implicitly, as already employed by Bjorken. 

At this point it should be noted that previous calcula- 
tions leading to results of the form (5) are, strictly speak- 
ing, inconsistent. On the one hand, an IR cut-off, which 
arises only from loop corrections to the tree level approx- 
imation, is needed in the calculation to screen the long- 
range interaction in the medium. On the other hand side, 
the (rcnormalizcd) vacuum parts of the same loop correc- 
tions (which make the coupling run) were omitted. This 
circumstance is the heart of the ambiguity that one has 
to guess a value of the parameter a in Eq. (5) - whereas 
a consistent calculation will actually specify the value of 
the coupling. 

Taking into account the running coupling (6) alters 
completely the structure of the integral (4), 1 
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Opposed to the former expression (4), the weighted cross 
section is manifestly UV finite - due to the asymptotic 
weakening of the strong interaction, large-i contributions 
are suppressed. For very hard jets, the integral (7) be- 
comes independent of the energy E, i.e., it is then con- 
trolled solely by the coupling at the screening scale. The 
necessary condition |ti| ~ TE 3> n 2 ~ m 2 D ~ aT 2 , is, 
for weak coupling, actually much less restrictive than the 
previous assumption £ > T to simplify the kinematics. 
Aspiring an extrapolation of the final result for dE/dx to 
larger coupling, where the hierarchy of scales \J~01T <C T 
breaks down, it is mentioned that the Debye mass is con- 
strained by tod Si 3T (see, e. g., [12]). Hence from this per- 
spective, the collisional loss could become ^-independent 
already for jet energies exceeding several GeV's - pro- 
vided, of course, that the present framework gives at least 
a semi-quantitative guidance at larger coupling, which will 
be advocated below. 

As an aside, it is instructive to illustrate in more detail 
the relation of Eqs. (4) and (7). The latter (full) result has 
the structure 
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which might be written in a 'factorized' form, 



In b — In a 
ln 2 r* 



ln 2 r* 



dr — . 

T 



where In t* = In a In b. With the notation a = 7 6, one 
can, if applicable, approximate In 2 r* = (In 6+2 In 7) In b w 
(In b + ln7) 2 , hence t* w 76 = \fab. Identifying b ~ 



1 Aiming here only at a leading order result, details related 
to the broken Lorentz invariance in the presence of a medium 
can be omitted. 
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TE/A 2 and a ~ [J? /A 2 , one observes that Bjorken's 'fac- 
torizcd' formula can be warrant - provided a is replaced 
by the running coupling a(t) at an intermediate scale 
ITE (which then leads to a cancella- 



t* ~ \ /'- / /■: ~ // / x 
tion with the logarithm in Eq. (4)). 

Returning to the full expression (7), which is actually 
both simpler and physically more intuitive, the collisional 
energy loss of a parton approaches 



dx 



q(m 2 ) 

b 



Ps(k) 



k 3 



(8) 



in the infinite-energy limit. With the ideal gas densities, 
the remaining integration yields a simple formula, 



dE n 



dx 
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(9) 



which differs clearly from previous expressions as Bjorken's 
formula (5). Aside from the modified cut-off dependence 
discussed above, the QCD collisional loss is proportional 
to a (instead of a 2 ). It is emphasized that the present 
considerations also show that the relevant scale for the 
coupling is the (parametrically soft) screening mass rather 
than a 'characteristic' thermal (hard) scale ~ T, as com- 
monly presumed. 

In order to quantitatively compare the new result to 
previous estimates it is necessary to specify parameters, 
namely A in Eq. (6) and the cut-off /i, i. e., the Debye mass. 
Similar renormalization arguments as above led in [12] to 
the conclusion that the familiar leading-order expression 
m 2 D = (l + AiraT 2 should be replaced by an implicit 
equation; to leading order 



= (1 



\ n f) 



A-Ka{m 2 D )T 2 



(10) 



whose solution can be given in terms of Lambert's func- 
tion. Obviously, this improved perturbative formula is jus- 
tified strictly only at temperatures T 3> A. Notwithstand- 
ing this, it is found in quantitative agreement with non- 
perturbative lattice QCD calculations down to T w l.2T c . 2 
I underline that the only 'free' parameter, fitted as A = 
0.2 GeV, is not only right in the expected ballpark (refut- 
ing the possibility of a mere 'coincidental' fit), but it also 
describes further non-perturbative results. As a matter of 
fact, the distance dependence of the strong coupling, as 
calculated from the heavy quark potential at T = in 
lattice QCD, is also reproduced quantitatively by the ele- 
mentary 1-loop formula (6) with t — > 1/r 2 - most notably 
even for large distances corresponding to a(r) « 1 [12]. 

It might be unexpected that resummed perturbation 
theory works beyond its expected range of validity 3 , but 
it is actually less surprising to find it appropriate jointly 

2 In contrast, the 'conventional' expression, again with the 
common choice a — > ot(2nT), deviates from the lattice results 
by an almost constant factor of 1.5 in a large temperature 
range, which is not much improved by the next-to-leading order 
correction. 

3 As argued e.g. in [13], perturbative expansions could (at 
best) be asymptotic series, implying an optimal truncation or- 
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Fig. 1. Light quark collisional energy loss: Eq. (9) vs. the 
prevalent expression (5) (which yields negative values for very 
large T) for representative jet energies. For details see text. 



for both quantities. In fact, the vacuum potential and 
screening at T > are closely related: they are both de- 
scribed by the (renormalized) i-channel scattering ampli- 
tude . . .which also determines the collisional energy loss. 
Put differently: one can renormalize the theory at T = 
(i. e., determine once and for all A from V qq {r)) 1 verify the 
applicability of the resummed perturbative approach for 
larger couplings as relevant near T c by successfully calcu- 
lating mrj(T), and give then justified estimates for dE/dx. 

With this motivation, an extrapolation of Eq. (9) as 
shown in Fig. 1 might be not too unreasonable. Assumed 
here is T c = 160 MeV and /i 2 = [\...2\m 2 D to estimate the 
uncertainty due to the IR cut-off. Shown for comparison 
are results from (5); here /i = mr) (likewise from (10)) 
albeit a in the prefactor (unjustified, but as commonly 
presumed) fixed at the scale Qt — 2ttT. It turns out that 
already near T c , the estimates from formula (9) exceed 
those from (5), even for rather large values of E. 

With regard to phenomenology it is instructive to take 
into consideration a principal effect of the strong inter- 
action near T c . From the distinct decrease of the QGP 
entropy established in lattice QCD calculations [14], one 
can, on rather general grounds, infer a similar behavior for 
the number densities p s . In the framework of the quasi- 
particle model [15], the particle densities in Eq. (8) read 
d s n±(y/ m 2 (T) + k 2 ), where the effective coupling in the 
quasiparticle masses, m 2 oc a e s(T)T 2 , is adjusted to lat- 
tice results for the entropy. The effect on the energy loss, 
as exemplified in Fig. 2, is sizable: near T c - despite the 
strong interaction - the QGP could become transparent 
due to the reduced number of 'active' degrees of freedom. 



der inversely proportional to the coupling. Hence for strong 
coupling, a leading order result might indeed be more ade- 
quate than sophisticated high-order calculations - somewhat 
opposed to naive intuition. 
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Fig. 2. The diminished scatterer density near T c reduces the 
energy loss, as shown here (for the case of gluons) by the dash- 
dotted curve. Below T » 1.2T C , Eq. (10) over-estimates the 
actual Debye mass (as obtained in lattice QCD); hence dE/dx 
could be slightly larger than depicted by the dotted line. 



It should be noted, however, that the detailed behavior 
of the (collisional) energy loss near T c is speculative. For 
instance, a distinctly dropping Debye mass, as suggested 
by some lattice calculations, could (partly) counterbalance 
the diminishing effect of the dropping particle density. 



4 Conclusions 

Reverting, for a resume, to a regime where perturbation 
theory can be expected viable, it has been demonstrated 
that taking into account the running of the strong cou- 
pling can influence the structure of results. For the colli- 
sional energy loss of hard partons, the result changes by 
one order in the coupling compared to previous expres- 
sions derived with a fixed a, cf. Eq. (9) vs. (5). While the 
increasing coupling at soft momenta leads to the para- 
metric enhancement, asymptotic freedom, on the other 
hand, implies that the collisional energy loss becomes in- 
dependent of the jet energy for E — > oo. Consequently, the 
asymptotic behavior of the underlying interaction makes 
the energy loss qualitatively different in QCD and QED 
(where an analog of Eq. (5) indeed holds). 

Except very near T c , Eq. (9) suggests a larger colli- 
sional energy loss than previously estimated [2,3]. This 
finding are a facet of the 'strongly coupled' QGP, which is 
characterized by large interaction rates. In fact, cross sec- 
tions a = f dtda/dt with running coupling are actually 
almost an order of magnitude larger than expected from 
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the widely used expression er Q fi x 



t)/V 



lmb. 



Thus, the present approach gives also a consistent expla- 
nation of the phenomenologically inferred large partonic 
cross sections O(10mb) [16,17] in the 'sQGP'. 



